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Abstract. The Hull-White one factor model is used to price interest rate options. The pa- 
rameters of the model are often calibrated to simple liquid instruments, in particular European 
swaptions. It is therefore very important to have very efficient pricing formula for simple in- 
struments. Such a formula is proposed here for European swaption. Based on a very efficient 
corrector type approximation the approximation is efficient both in term of precision and in term 
of spped. In our implementation the approximation is more than ten time faster than the direct 
pricing formula and more than twenty time faster than the Jamshidian trick. 



1. Introduction 



Due to its flexibility, the iHull and White! (|1990f ) one factor model is often use as a simple and 



efficient model to price interest rate derivatives, including exotics. 

The model parameters are often obtained through a calibration procedure. A calibration is 
a reverse engineering work were the model parameters are reconstructed from market prices. In 
most of the cases this inversion is not obtained explicitly but implicitly by solving numerically a 
non-linear equation. The price function being the one to inverse. In the numerical procedure, the 
pricing function is used repeatedly. This is the reason why a numerically efficient price computation 
of vanilla instruments is very important. 

In this note the price of European swaption in the Hull- White one factor model is dis cussed. The 



first e xact pricing solution proposed for that model is probably the one proposed bv I Jamshidian 



(1989). Its solution is based on a decomposition, now called Jamshidian's trick. The decomposi- 



tion consists in dividing the bond (or swap) in a set of zero-coupon bonds with strikes such that 
all of them arc- exercised in the same conditions. As the model is a one-factor model with prices 
monotonous in the underlying stochastic variable, such a level exists. In the Jamshidian decompo- 
sition, the underlying random variable is the short-rate. A one dimensional equation consisting of 
exponentials with relatively complex coefficients needs to be solved to find the (common) exercise 
level. The pricing is obtained by adding zero-coupon options. For a n coupon bond, there are 2*n 
normal cumulative functions to compute. 



A more direct solution was proposed in iHenrardl (|2003l ). Using an explicit description of the 



zero-coupon price in term of a normally distributed random variable, the swaption price is written 
as the discounted cash-flow multiplied by some weights (probabilities). The exercise level in term 
of the random variable is computed by solving a one dimensional equation. The equation is similar 
to the one in the Jamshidian decomposition; it consists of exponentials but its coefficients are more 
easy to write. The price is obtained with one normal cumulative functions for each coupon plus 
one for the strike. For a n coupon bond, there are n + 1 such functions to compute. Depending of 
the implementation, this second method is usually two to four times faster than the Jamshidian 
method. 
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2. Model, hypothesis and preliminary results 



The HJM framework describes the behavior of P(t,u), the price in t of the zero-coupon bond 
paying 1 in u (0 < t, u < T). When the discount curve P(t, .) is absolutely continuous, which 
is something that is always the case in practice as the curve is constructed by some kind of 
interpolation, there exists fit, u) such that 



(1) 



P(t, u) = exp 



f(t,s)ds 



The idea of Heath et al. ( 1992h was to exploit this property by modeling / with a stochastic 
differential equation 

df(t, u) = fx(t, u)dt + a(t, u)dW t 

for some suitable (potentially stochastic) fi and a and deducing the behavior of P from there. To 
ensure the arbitrage-free property of the model, a relationship between the drift and the volatility 
is required. The model technical details can be found in t he original paper or in the chapter 
Dynamical term structure model of lliunt and Kennedy ( 2004 ). 

The probability space is (O, {J~t}, T, P). The filtration Tt is the (augmented) filtration of a 
one-dimensional standard Brownian motion (Wt)o<t<T- To simplify the writing in the rest of the 
paper, the notation 



vit, u) 



a(t, s)ds 



is used. 

Let Nt = exp( J Q r s ds) be the cash-account numeraire with (r s )o< s <T the short rate given by 
r t = f(t,t). The equations of the model in the numeraire measure associated to Nt are 

df(t, u) = a(t, u)u(t, u)dt + a(t, u)dW t 

or 

dP N (t,u) = -P N (t,u)v(t,u)dW t 

The notation P N (t, s) designates the numeraire rebased value of P, i.e. P N (t, s) = iV t " 1 P(i 3 s). 
The following separability hypothesis will be used: 

H: The function a satisfies a(t, u) = g{t)h{u) for some positive functions g and h. 

Note that this condition is essentially equivalent to the condition (H2) of iHenrardl ( 2003f ) but 
written on a instead of on v. The conditi on on v was v{s,t2) —jAs,t\) = f(ti,t2)g(s). 
We recall the generic pricing theorem ( Hunt and Kennedy . 2004L Theorem 7.33-7.34). 



Theorem 1. Let Vr be some Tt -measurable random variable. IfVr is attainable, then the time-t 
value of the derivative is given by V t N = V N + J * <f> 8 dP? where fa is the strategy and 

V t =N t E N [V T N^ l \T t ] . 

We now state two technical lemmas that were presented in IHenrardl (120051 ). 



Lemma 1. Let < t < u < v. In a HJM one factor model, the price of the zero coupon bond can 
be written has, 



P(u,v) 



P(t,v) 



exp 



f 

'2.,, 



P(t,u 

Lemma 2. In the HJM one factor model, we have 



(y (s,v) — v 2 {s,uj) ds + / (i/(s,v) — u(s,u)) dWi 



N U N~ 



exp 



(- J r s ds^j =P(u,v)exp^y v(s,v)dW s - ^ J v 2 {s,v)ds 
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We give the pricing formula for swaptions for a future time (see ( Henrardl . 20061 . Theorem 2)) 



The notations used to describe the swaption are the following. Let 9 be the expiry date and the 
swap is represented by its cash-flow equivalent (ti, Ci)i=o,..., n - The date to is the swap start date 
and U (i = 1, . . . , n) are the fix coupon dates. The amounts Co is —10, c% > (i = 1, . . . , n— 1) are 
the coupons and c n > is the final coupon plus 1 for the notional. 

The swap is described by its cash-flow equivalent (ci,ti) (0 < i < n). The dates ti are the 
settlement date to, the (fixed) coupon payment dates ti (1 < i < n). In the descriptions a receiver 
swaption is considered. The values are the initial notional co = — 1 (resp. strike Co = —K), the 
different coupons Ci and the final cash-flow including the coupon and the notional. Usually the 
coupons are regular with Ci = SiRk for a strike rate Rk and a day count fraction Si between ti 
and ij+i. The swaption expiry is denoted 9 and 9 < to. 

Theorem 2. Suppose we work in the HJM one-factor model with a volatility term of the form 
(H2). Let 9 < to < • • ■ < t n , cq < and Ci > (1 < i < n). The price of an European receiver 
swaption, with expiry 9 on a swap with cash-flows Ci and cash-flow dates U is given at time t by 
the Tt-measurable random variable 

n 

^QP^t^iV^ + 

where K is the Tt-measurable random variable defined as the (unique) solution of 



(2) 5jciP(t,ti)exp -- 



-a? - ctiK I = 



i=0 

and 

2 _ / /,.(„ + \ „r„ / n^^2 



a/ = J {v(s,ti)-v{s,9)yds. 
The price of the payer swaption is 

n 

-^ / C l P(t, ti)N(-K-Oi) 
i=0 

The Brownian motion change between the N t and the P(t, tj+i) numeraires is given by 

dW 3 t +1 = dW t + v{t, t J+1 )dt. 

3. Swaption pricing 

From an option pricing point of view, swaptions and bond options are equivalent. From now 
on we will refer only to swaption. Up to the change of names, all the results are valid both for 
swaptions and bond's options. 

The forward values of the zero-coupon bonds and swap without the initial notional are given in 

and Bt = _%__%M. 

* P(Mo) P(Mo) 

Those value are the value rebased by the numeraire P(., to) 

Let v % (t) = (y(t,ti) — v(t,to)) (0 < i < n). In the martingale probability associated to the 
numeraire P(.,to), the rebased prices are martingale that satisfy the equations 

dP^-Py(t)-dW t °. 

The zero-coupon bond prices are exactly log-normal as the volatility v l (t) is deterministic. 



4t is — K for a bond option of strike K. 
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The rebased swap value satisfy 



dB t = - J2 c t p iv l (t) ■ dW°. 



i=l 

Using the notation a\ = CiPl/B t and a(t) = ^ a\v l (t), the equation is becoming 

dB t = -B t a{t)- dW°. 

This is formally a log- normal equation but the a coefficient is state dependent. 
The volatility a(t) of the coupon bond can be approximated by its initial value 

P 



Let 



Mt)|d*. 



With the approximation the equation for the bond B is log-normal and the standard Black formula 
approach can be used: 



Ro = P(0,t o )E° p- 1 (0,t o )(S e -ifF(0,to)) + ^P(0,t o )E° 



<7 X - i(7 Z 



Theorem 3 (Initial freeze approximation). In the HJM model, the price, with initial freeze ap- 
proximation, of a receiver swaption with expiry 8 is given at time by 



Ro = P(0, t ) {B q N(k + <f ) - KN(k)) 



where 



1 

Til 



]n(B /K)-^a ' 



The price of a payer swaption is 

P = P(0, t ) (KN(-k) - B q N(-k 



*o)) 



This formula is equivalent to pricing the bond with the Black formula and an implied volatility 
eq ual to dn . 

Baviera (120061) proved that this approximation is very efficient in the context of the BMM. 
Intuitively the largest part of the swap is in the notional at maturity t n . That part is exactly 
log-normal and the total is well approximated by a log-normal dynamic. In the example analyzed 
in the relevant section it will showed that this is true also for amortized swaptions for which the 
notional is paid by in trenches along the swap's life. 



3.1. Corrector style approximation 

Carlo simulations linked to the LMM. It was initially developed by 



The predictor-corrector method is often used in M onte 

Kloden and Platen! (| 19951 ) and 



used in the framework of the Libor Market Model bv lHunter et al 



(|2001l ). It consists for a given 



path to first simulate a step with the Euler scheme (equivalent to the initial freeze approximation) . 
The value obtained is used as an approximation of the step end state (predictor). Then the (state 
dependent) coefficient is approximated by the average of the coefficient at the starting state and 
the coefficient at the approximated end state (corrector). This approach leads to very efficient 
Monte Carlo implementations. 

The above described method is implemented at the path level. For each path the approximated 
final value is different and the state dependent coefficients approximation is different. The path 
level approximation is suitable for Monte Carlo implementation but not for approximated explicit 
formulas. It is possible to design a method that works with the same idea of averaging coefficients, 
not at the path level but at a more general level. The standard approach works at the path level 
and is the same for all options priced with the same Monte Carlo. The new approach described 
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here extends the one proposed in iHenrardl ( 2007 ) and is based on approximation adapted to the 
option priced and valid for all paths. 

The approach is first described in an intuitive and loose way before going through the details. In 
the Black- Scholes formula the option price can be viewed as today's prices of the option components 
(Bo, K) weighted by some probabilities. If the probabilities are well approximated, the price is 
good. The main probability used is the one of reaching the strike from the current position. In 
the Black-Scholes formula it is roughly N(d\). For the paths that reach the strike the coefficients 
can be approximated in the corrector way. The corrector state is given by a strike state. The 
corrected volatility is the average between the initial coefficient and the strike (final) one. In our 
case the figure represents the coefficient of a log-normal distribution that ends up on the strike for 
the same paths as the real distribution. Said in another way, this is the (approximative) equivalent 
log-normal volatility for which the probability to reach the strike is equal to the actual one. In 
options term, it is the implied volatility that gives the same price for that strike. 

The value of the different rates and zero-coupon bond prices at the strike are not unique as 
the model is multi-factors. The exact way the strike information is transformed into a coefficient 
information need to be decided. There is one constraint for m parameters. The choice is within 
a 77z — 1-manifold into a m dimensional space: an infinity of potentially large dimension but with 
zero-measure. 

The description of the option price in term of probability is not exact. The option is a probability 
average of the pay-off. The price obtained through the approach described above will not be exact. 
It is a significant improvement to the initial freeze/Euler approach. 

A strike value of the different parameters has to be selected. The swap price is at-the-money at 
expiry when 



E< 



i=0 

The discounting value of the zero-coupon bond can be approximated (initial freeze) by 

9 P(0,£ o ) V 2 



with 



rf= / y (t)\ 2 ds 
Jo 



and the Xi standard normally distributed random variables. By choosing (arbitrarily) to have all 
the stochastic variables Xi equal at the strike the (one dimensional) equation to solve is 



1 2 



=0 



^ c t P£ exp [ -nx - -Tj ) = 0. 



Obtaining the solution to the above equation requires to solve a one dimensi onal equat i on equ ivalent 
to the one solved in the swaption price for the Gaussian HJM model (jHenrardl (|2003f )). For 



numerical reasons one may prefers not to have to solve this type of equation. The above equation 
can be replaced by it first order approximation 



E 

i=0 



,.Pq I - 



2 



the solution of which is explicit: 



V™ r P l - IV™ r-P i T 



The zero-coupon prices, in the exponential case and the approximated first order case, are given 

by 

p k = K cx P ( ~nx - o r n i respectively P l K = P^ ( 1 - nx - -t? 
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The rates and bond prices are 



Y[(l + 5 ] Li K ) = (Pk)- 1 and B K = Y,c i P A K = K. 

3=0 i=l 



By defining 



1 K 



CiP % K 



B K 

the swaption can be priced with a (option strike dependent) approximated volatility 
(3) VK{t) = \(^{a + tf K )v\t)\. 



Note that in the multi- factor model, the volatility (JK{t) is a vector, as vq and vk are. 

Theorem 4 (Corrector implied volatility). In the local volatility LMM, the price, with initial freeze 
and corrector approximation, of a receiver swaption is given at time by 

Ro = P(0, t ) (B N(k k + <J K ) - KN(k k )) 

where 

K K = — (ln(B /K)-la 2 K 

0~K \ 2 



°k= / \v K {t)\ 2 dt. 
Jo 



The price of a payer swaption is 

P = P(0, t ) {KN{-k k ) - B N(-k k - a K )) 

3.2. Hull-White volatility. The different coefficients used in the previous sections are explicitly 
described for the Hull- White one factor model. 

In the case of the extended Vasicek or one-factor Hull and White model, one has a(s,i) = 
7yexp(— a(t — s)) and v(s,t) = (1 — cxp(— a(t — s)))r]/a in the constant volatility case. The time- 
dependent volatility is also covered with a(s, t) = rj(s) exp(— a(t — s)) and v (s, t) = (1 — exp(— a(t — 
s)))rj(s) / a. The a used in Theorem [2] are given in the constant volatility case by 

«? = ITS (exp(-afl) - cxp(-at J )) 2 (exp(2a(9) - 1)) . 
2a J 

In the time-dependent case, rj is piece- wise constant with r/(s) = rjt for Sj_i < s < s,; and 
= so < Si < • • • < s n = +oo. The expiry dates are between some of those dates and the 
relevant dates are denoted s g _i < < s q . To shorten the notation an intermediary notation is 
used: ro = sq < r; = s; < r q = 9. With those notations, one has 

9-1 



«? = ttt (exp(-a0) - exp(-aii)) V" r;, 2 (exp(2ar;+i) - exp(2an)) . 
The zero-coupon bond instantaneous volatility v % (t) is given by 



v l (t) — — 7](s) exp(ai)(exp(— ato) — exp(— atj)). 
a 

The [0, #]-period volatility is similar to ai and is given by 

9-1 



1 I > 

— (exp(-aio) - cxp(-a^)) 2 ^ Vi (exp(2an+i) - cxp(2ar ; )) 



^2 _ 

~ 2a 3 



1=0 

From there all the coefficients used in the pricing follow from simple arithmetic operations. 
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4. Approximation numerical analysis and comparison 

4.1. Approximation quality. The approximation quality is estimated in the following way. An 
extended Vasicek model with time volatilities is calibrated to market data. The data are the 
yield curve and the at-the-money volatilities. This is done with two different date and to types of 
calibration to ensure some diversity in the test. 

Once the model is calibrated, it is used to price a wide range of swaptions using both the exact 
pricing formula and the approximated formula. The Black swaption implied volatility of the results 
is used as comparison figure. In term of precision, the market quote the Black volatilities with a 
maximum precision of 0.10%. Any approximation error below this can be consider as below the 
measure error. Nevertheless we a priori require a error below 0.025% for options up to 300bps 
out-of-the money. It is compared to other approximations proposed in the literature. 

The first example is done with market data from 4 July 2008. The model is calibrated to 
diagonal ATM options with expiry between 1 and 10 years and tenors between 10 and 1 year. The 
mean reversion parameter used is 0.02. 

With that calibration, four type of swaptions are priced: 1Y x 10Y, 5Y x 5Y, 8Y x 2Y and 2Y x 
20Y. For each of those types, a large range of strikes, up to 300 bps away from the money is used. 
The results are given in Figure [T] 




-1 o 1 

Strike vs ATM 

(a) Volatility 




Strike vs ATM 

(b) Price 



Figure 1. Error introduced by the approximation. Market data as of 4 July 2008. 



The Figure contains the error for the proposed approximation and the initial freeze approxi- 
mation. The improvement of the method is clear on the picture. For all the tested swaps, except 
for the longest maturity and the most extreme strikes, the error is within the a priori band. Note 
that the options for which the error is above the a priori limit the strike is 300 bps away from the 
money; that level of moncyness is usually not quoted in the market. For most of the errors, it 
is difficult to distinguish them from in the graph. For the most extreme case, where the error 
is above the a priori limit, it may be useful to look at the error in term of price instead than in 
term of volatility. For deep out-of-the-money option an error on the volatility will have a smaller 



impact on the price (Figure 1(b) I. In that case it appear clearly that even if the extreme strikes 



have higher volatility error, the price error is relatively small. 

The first example showed that the error is larger for long tenor options. For that reason the 
second example is concentrated on those options. The results are reproduced in Figure [2] Note 
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that in the second example, the strike below the money are only up to 250 bps away from the 
money. As the rates are around 3%, a 300bps shift would give a strike close to 0. 

In Figure 1(a) the proposed approximation is compared to other app roximation prop osed in 
the literature. The first comparison is with the initial freeze proposed by iBavieral (j2006l i ( in the 
context of the BMM). In the graph, the initial freeze is represented by the dashed lines. That 
approximation is efficient for short tenor options but relatively inef ficient for long tenors (10, 20 
years) and strikes more than 100 bps from the money As described in lrlenrardl ( 20081 ) the proposed 
approximation is at least ten times more efficient. 

The second comparison i s with the low variance martingale (LVM) approximation proposed by 
Schrager and Pelsserl (|2006l ) . The approximation is represented in the graph by the dotted-dashed 
lines. The quality of that approximation is generally below the initial freeze one. The error is 
above the 0.10% limit as soon as one goes more than 100 bps away from the money. The figures 
reported here are in line with the one reported in Table 4.2 of the original paper. The error is 
generally more than 25 times larger that the one with the approximation proposed here. 





(a) Volatility 



(b) Price 



Figure 2. Error introduced by the approximation. Market data as of 19 decembre 2008. 

The second example, with market data significantly different, confirm the result of the first test. 
Up to 10 years tenors, the proposed approximation is sufficient. For longer tenors and strikes very 
far away from the money (more than 200 bps), the error may be above the a priori very stringent 
limit of 0.025% but still below the market precision of 0.10%. Note also that at that level of strikes, 
the market is relatively illiquid and their is a lot more uncertainly in the market level than in the 
approximation. 



4.2. Speed. The speed testing is done using the last swaption in the second test. The price is 
computed on the author computer using a Matlab implementation. The price is computed 1000 
times and the time required for the computations is recorded for each of the three implementations 
(exact, corrector approximation, initial freeze approximation). The result is 9.73s in the first case, 
0.62s in the second and 0.58 in the third. The corrector approximation is around 15 times faster 
than the exact formula. The difference is in the computation of the k, which is obtained by solving 
a one-dimensional non-linear equation. This takes around 90% of the time for the exact formula. 
The speed gain using the initial freeze is only 8% with respect to the corrector. Given the error 
which is significantly larger, the initial freeze approach is not advisable for calibration. The LVM 
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approximation takes around 50% more time than the corrector approximation and its precision is 
far less. There is no reason to use it in the Hull- White one-factor context. 

5. Conclusion 

An approximation for the pricing of European swaption in the Hull- White/extended Vasicek 
one factor model is proposed. The approximation is based on a corrector style approach. The 
approximation is proved very efficient with an error below the market precision. The approximation 
quality decreases with the underlying swap tenor and the swaption moneyness. Only for very 
long tenors and very far away from the money options the approximation quality decreases to be 
noticeable but always inferior to market precision. 

The time efficiency of the implementation with corrector approximation is largely superior to 
the one of the exact formula. In our implementation the speed is increased fifteen fold. When 
a large number of computations is required, like in calibration, this gain can become significant. 
When speed is required we suggest to use the corrector approximation implementation. 

The approximation approach presented in this note was initially developed for multi-factor 
Libor Market Model. It is more generic than presented here. It can also be applied to multi-factor 
Gaussian HJM model like the G2++ model. 

Acknowledgment: This note is dedicated to an efficient person who's birthdate is on December 
20th, the date this note was started. 

Disclaimer: The views expressed here are those of the author and not necessarily those of its 
employers. 
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